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Abstract 

The intrinsic length scales of a reaction diffusion system (Gierer-Meinhardt 
model)is varied by quasi-statically changing the diffusion constant of the activator 
and a transition from rolls to hexagon is detected. The transition is hysteretic or 
first order like. From stability analysis, we also analytically show the possibility of 
such transitions. 

PACS number (s): 87.10.+e, 47.70.Fw 

Multilateral structures like squares, hexagons, rhombs are observed in chemical pat- 
terns [1], in Faraday surface wave ([2, 3] and references therein), nonlinear optical systems 
[4] and in many other systems. Many reaction diffusion systems are studied to show 
such multilateral structures under various ambient conditions. The origin, stability and 
dynamics of such multilateral structures are subjects of active investigation [5, 6, 7, 8]. 
In a reaction diffusion system the diffusivities are taken to be constants of activators and 
inhibitors and the variations in them due to variations in temperature or other experimen- 
tal conditions are considered to be on very small scales to produce any effect. We argue 
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that if there are some sort of phase transitions which arc hystcrctic with the variation of 
diffusion constant, then it indeed has relevance to study transitions under the variation 
of diffusivities. In a situation when diffusion constants are very close to such transition 
points, a slight change in them due to the fluctuations in temperature etc. of the system 
might bring about a phase transition which will not be reversed when the diffusivities are 
restored. In another respect such studies are worth doing is to understand how such a 
system behaves under a slow variation of its internal length scales. A slow variation of 
intrinsic length scale can be effected by slowly varying the diffusion constant of the sys- 
tem. In the present work we are going to represent and characterize the phase transitions 
brought about by the variation of diffusivity of the activator of Gicrer-Meinhardt (GM) 
[10] model keeping the diffusivity of the inhibitor fixed. In what follows, we will show roll 
patterns are loosing their stability to hexagonal structures via a first order like transition. 
Under the action of competiting length scales, the closed multilateral structures charac- 
terized by wave numbers which locally satisfy the condition qi + q2 + + = are 

indeed an impressive general outcome and regular hexagons are a more symmetric special 
case of them. 

The Gierer-Meinhardt model that we have taken up is 



Where A is the activator concentration and B is that of the inhibitor species. In what 
follows, we will always take the Turing [9] condition Db > between the diffusivities 
of activator and inhibitor to hold true. In the above expression pA ^ Pb are the reaction 
strengths, ha & A*b are the self removal rates and & are the basic production 
terms for A and B respectively. In a parameter region where steady spatial order of 
concentration forms the relaxant length scale is given by [11] 



where Ci and C2 are constants. Thus the diffusivities are the one which selects the internal 
length scale of the system. The present report comprises of three parts. The next part 
shows the results of numerical simulation done under quasi-static variation of the diffusion 
constant Da- The following part accounts for the observations of numerical simulation 
on the basis of a stability analysis. The last part includes discussion. 

In our simulation we have set aA — ctb — 0. The other parameters are pA = Pa = 0.01, 
Pb — l^B — 0.02, Ka — .25, Db — 0.2 and Da is varied from to 0.01 by the steps of 
0.00001 near the transition regions and by the steps of 0.0001 in other parts keeping 10^ 
temporal steps between every variation of Da- The simulation is done on two dimensions 
by implementing finite difference method and on a lattices of size 256 x 256. In the first 
part of our two dimensional simulation so long as Da remains less than approx. 0.004, 
the locally isolated concentration peaks show up. Such structures appear for very small 
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value of Da compared to Db because of the fact that the activator finds no time to spread 
compared to very fast inhibitor spreading. As Da is increased, near Da = 0.004 steady 
roll structures which are locally parallel appear and persist upto Da < 0.008. As Fig. 1(a) 
shows, many small domains of parallel rolls are present and are oriented in all possible 
directions. These domains are separated by very many defects and domain walls. At 
about Da = 0.008 a transition from rolls to almost regular looking hexagonal structures 
takes place. This phase is stable upto a Da value about 0.009. On a further increase 
in Da, hexagonal patterns loose its stability to a steady homogeneous state. Fig. 1(a) 
and 1(b) show the roll structure at Da — 0.005 while increasing Da and hexagonal at 
Da — 0.0085 as seen in the lattice. After the hexagons have stabilized at Da > 0.008, 
now if we go on decreasing Da the rolls do not come back at Da = 0.008 but at a much 
lower value. The transition from rolls to hexagons show hysteresis. A plot of average 
concentration of A against Da is shown in Fig. 2. These plots clearly show a jump in 
the average concentration as we go from below (circles) along Da axis which marks the 
transition to hexagons from rolls at Da = 0.008. As we come down the Da axis (triangles), 
we see that the average amplitude does not immediately jump back to the lower line at 
Da = 0.008 but continues to follow the same upper line representing the hexagons and 
comes back to the lower value at around Da = .0055 which marks transition to rolls from 
hexagons. Thus Fig. 2 shows hysteresis and identifies the transition to be first order like. 
The average activator amplitude also comes out as a relevant parameter to characterize 
such situations. 

Having numerically found the results let us try to explain it from analytics. After 
rescaling of concentration, length and time as A = /jaA, I = \J and t = iixt respec- 
tively, Eq.l takes the form 

: D^^A + pA- 
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where we have taken ua — cfb — 0. Under such scaling, D — ^ < 1, pA — 

K — Pb — and Ub — The homogeneous steady fixed point of the system is 

Ma Ma Ma 

* 2 

* * - A 

A —A— .0085 (under prevalent conditions of the simulation) and B =B= — 85.0. A 
linear stability analysis shows that inhomogeneous perturbations will grow with a growth 
rate A given by 
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In the range of given parameter values the real part is always negative giving no scope 
for Hopf modes to appear and in our simulation we have not seen any too. For a roll 
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solution to grow the discriminant has to be greater than zero and its absohite vahie has to 
be more than the first part in the expression of A. These two conditions give the following 
inequalities for the range of value of the selected wave number with the variation of Da 
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These graphs are plotted on a A; — D plane as shown in Fig. 3. In this figure we can 
easily identify the region where rolls will grow as that bounded by the two curves and the 
k axis. In the above expression k is the magnitude of k. In our simulation this region has 
been captured for Da > 0.004 to 0.008. 

The interesting things show up if we try to see how a perturbation grows which tend 
to rotate the existing roll solutions. To do so we perturb the roll solutions as A{1 + 
(5acos (ki.r)) cos (k.r) so as to probe a possible rotation of the locally parallel rolls and 
we end up with two sets of linear decoupled equations in 6a and db obtained from harmonic 
balance of the terms containing cos (ki + k).r and cos (ki — k).r. The equations are 
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It is apparent from the shape of above two sets of equations that a simultaneous 
vahdity of them requires ki to have two values such as ki = kn and ki = ki2 where 
|kii I = |ki2| and kn & ki2 are at angles 9 and n±9 with the direction of k. Such a choice 
of ki predicts development of rhombic structures at the cost of rolls. But we have not 
been able to stabilize such structures. A regular hexagon grows when |kii| = |ki2| = |k| 
along with |ki| =0 and 6 = 60. Because in such a situation we get all the three wave 
numbers, one being that of the preexisting roll whose growth rate is the same as other 
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two new modes which are locally rotated at 60 with it, grow equally in the concentration 
of the activator A obeying the equation as follows. 
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. Now putting the relative values of the 



eigen values corresponding to A and B we get an inequality specifying the wave number 
arne for the growth of hexagons as 
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This curve is also plotted on Fig.3 where it almost superposes on the graph given by the 
Eq.(5) and hexagonal order should show up in the region below this curve. Thus the 
hexagons and rolls can coexist over a considerable region of this phase space. In fact in 
almost all the regions where rolls are stable hexagons are also stable. This fact might 
account for the fact of obtaining so many small domains of locally parallel rolls and the 
abundance of defects. As we go on increasing Da, we eventually come to a region where 
the Eq.(4) is violated and thus rolls become unstable to hexagons for which we are still 
in the stable region. It is important to note that Fig. 1(b) shows a field of hexagons 
with very few defects compared to those in Fig. 1(a). This is because the hexagons are 
formed in a region where rolls are not stable. A further increase in takes us beyond the 
boundary given by Eq.(8) and hexagons finally loose stability to a homogeneous steady 
state. With this phase diagram we cannot directly account for hysteresis at the transition 
region of rolls and hexagons but a wide region of overlap is quite indicative of the fact 
that such things can happen. 

In the conclusion, we would like to mention a transition from rolls to hexagons has 
been observed with increasing diffusivity of the activator species in the GM model. The 
transition region shows hysteresis. Such hysteresis can cause a permanent change in the 
phase of the system even with a small fluctuation in the diffusivities which are generally 
treated as constants. We also analytically captured the fact that closed multilateral 
structures can not simultaneously show up in the concentrations of activator and inhibitor 
and the system decouples in such a situation. At a very high diffusivity of the activators 
the hexagons looses stabihty to homogeneous steady state. We argue that a breakdown 
of Turing condition might be the cause of such a phase change. 
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FIGURE CAPTIONS 

Figure 1(a).- The figure shows many domains of locally parallel rolls separated by domain 
walls and other defects at a value of the diffusivity Da — 0.005. In this figure the 
concentration of the activator (A) has been plotted on a 256x256 lattice space. 



Figure 1(b).- The figure shows regular hexagonal order in the concentration of the 
activator species A at a value of the diffusivity Da — 0.0085 on the same space lattice as 
Fig. 1(a). 



Figure 2.- The figure shows a plot of the average activator concentration against Da 
where the circles arc plotted for increasing Da and the triangles arc plotted when Da is 
decreasing. At Da = 0.008, a jump in the curve indicates transition to hexagons from 
rolls and at Da = 0.0055, the fall in the curve marks coming back from hexagons to rolls. 



Figure 3.- The figure shows the phase diagram as obtained from linear stability on a 
D — k plane. The roll pattern is stable in the region above the curve '1' and below the 
curve '2'. The hexagons are stable anywhere below the curve '3'. Thus there is a region 
of coexistence for rolls and hexagons which he above curve '2' and below curve'3'. In the 
region below curve '2' and at the same time below curve '3' only hexagons are stable. 
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